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Abstract 
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results to theoretical predictions. 
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I. INTRODUCTION 



Quantum phase transitions (QPTs) describe sudden changes of the ground state of a 
many-body quantum system as a non-thermal control parameter moves through some criti- 
cal value l| (at zero temperature). QPTs are relevant not only for understanding of quan- 
tum many-body systems, but also for other problems such as quantum entanglement ^ 
and quantum computing, e.g., adiabatic quantum computing 3[ and quantum estimation 
4j. Interesting phenomena related to QPTs have recently been experimentally observed in 
various systems, such as heavy fermions and Bose-Einstein condensates j^]. 

There has been a recent flurry of activity following the observation jf| that the proximity 
to a quantum critical point enhances the sensitivity of a system to external perturbations, as 



measured by quantum-information-theoretical quantities such as the Loschmidt echo 



6j or 



the ground state fidelity 3]. Exploiting such sensitivity, one can detect quantum criticality 
by coupling an additional spin as a probe to the system undergoing a QPT. This was 
suggested in and demonstrated in ^j], where the local coupling to the probe qubit was 
used as the perturbation. 

Here, we implement an alternative method to detect the critical point of a QPT by 
measuring an arbitrary qubit of the quantum critical system while applying a global per- 
turbation. The critical parameters of a general QPT, i.e., including critical points and 
exponents, can in principle be detected using our method. Our approach does not require 
an additional probe spin, which makes the experimental implementation easier. In contrast 
to our method, in the previous approach 9| the efficiency of detection depended on the 
nature of the phases on both sides of the critical points, and could be affected, or even 
rendered insensitive, by the locality of the probe. Furthermore, because our method uses a 
global perturbation, it increases the echo signal - making it, in principle, better suited for 
scalability with the size of the system. 

The paper is organized as follows: In section II we introduce the model and discuss 
how we use it to simulate a second order QPT. In section III we review the behavior of the 
Loschmidt echo in a critical system using a perturbative treatment. In particular, we discuss 
the echo decay rate and its scaling near the critical point. In section IV we describe the 
experimental implementation for even and odd spin chains using nuclear magnetic resonance, 
and compare our results to theoretical expectations. We offer concluding remarks in section 
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V. 



II. ISING CHAIN WITH A TILTED FIELD 



To demonstrate the detection of quantum criticality, we simulate the QPTs using a one- 
dimensional antiferromagnetic Ising model with the Hamiltonian 

N-l N N 

i=l i=l i=l 

where B z and B x denote longitudinal and transverse magnetic fields, respectively, <j\ and 
a l x are Pauli matrices acting on spin i of the chain, and the coupling strength has been set 
to unity. This type of model has been extensively studied in the literature in the contexts 



of statistical 
15 



16. 
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rysics [101 



ia 
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quantum computing [ll| . quantum chaos [12], and QPTs 



22|. 



HQ, 



Notice that the general case of Eq. ([!]) with B z ^ and B x ^ cannot be solved 
exactly using Jordan- Wigner transformation methods because the longitudinal field maps 
into high order coupling of the resulting fermions. This can also be seen by noting that the 

, with B z the longitudinal 



Hamiltonian ([T]) can be mapped into a classical 2D Ising model 
field and B x an effective temperature - which means that our quantum simulation can also be 
seen as a simulation of this archetypal model of classical phase transitions. The map between 



a quantum d dimensional spin system into a d+1 classical Ising system 



231 ] lets us obtain the 



phase diagram of Hamiltonian ([!]) in the thermodynamic limit, which corresponds to that 



of the 2D classical antiferromagnetic Ising model [24|, and is shown qualitatively in Figure 
[TJ The critical line is second order except for B x = 0, where it is a first order transition. As 
we will use only a few qubits, we are concerned here only with finite systems. Furthermore, 
the Loschmidt echo decay rate typically increases with system size [8j, which implies that 
in the thermodynamic limit the echo would decay infinitely fast (unless the perturbation 
is simultaneously reduced to zero, where a singular decay rate would be obtained Q). In 
the finite size systems under consideration, the gap across a second order transition never 
closes, but rather reaches a minimum near the critical point (this minimum goes to zero in 
the thermodynamic limit). Furthermore, for finite systems we need to consider odd-even 
effects, which in our model system will introduce "quasi" -phases that come from boundary 
effects and merge in the thermodynamic limit. 



Let us consider first the ground states for B x = 0, which will be relevant for our exper- 
iments. We keep in mind that in this particular case the system undergoes crossovers as a 
function of B z , since only the energies, not the eigenstates, depend on B z . When N is an 
odd integer, the ground state of the system is 



\V{B Z )) 
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0) (-2 < B z < 0) 



(2) 
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where |0) and |1) are the eigenstates of a z . We denote the four phases of the ground state 
as \ipl) with k — 1, . . ., 4. The energy of the ground state is 
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We denote the energy corresponding to the four phases \if)%) as E° k , respectively. B z = 
B c = ±2 and are the crossover points, where the system has a degenerate ground state. 
(N + 1 )/2 degenerate states exist at B c = ±2, making them the multiphase points of the 



system 



15 



25]. 



When tV is an even integer larger than 2, the ground state of the system is 
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and the energy of the ground state is 
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The crossover points are B c = ±2 and B c = ±1. Points _B C = ±2 are also multiphase points, 
each with iV/2 degenerate states. The five phases are denoted as with k — 1, . . ., 5 and 
the corresponding energy is represented as E e k . 

From Eqs. © and ©, one finds that if iV -> oo, £° 2 -> ££ 3 , ££ 2 -> ££ 3 , and 
E^ 4 — > -E 1 ^ 3 - Hence, in the thermodynamic limit only the multiphase points B c = ±2 are 
the crossover points, and l^), |^ 2 )' iV'f); an d l^f) are "quasi" -phases that merge into 
a single antiferromagnetic phase (see Figured]). The finite size energy phase diagrams are 
shown in Figure [2](a-b). 

In general, second order QPTs are characterized by a closing of the gap between the 
ground and first excited energy levels at the critical points (in the thermodynamic limit). 
Using our small quantum information processors, we will simulate the evolution of the 
quantum system described by Hamiltonian (CQ) in a regime where its spectrum is similar to 



the general case of a finite-size second order QPT (that is, with a small but finite gap). 
We achieve this by using a small transverse field B x to lift the degeneracy at points B c , 

n 

which makes the spectra resemble a continuous QPT [26(. Thus, we explore the transitions 
crossed by the dashed line in Figure [TJ In the analysis of our results we must take into 
consideration finite-size effects such as the size of the gap at the critical points, and the 
additional "quasi" -phases introduced by flipping a finite number of spins at the ends of the 
chain - which makes a distinction between experiments with odd and even chains. 



III. THE LOSCHMIDT ECHO AND QUANTUM PHASE TRANSITIONS 



A. Detection of critical parameters 

Let us consider a system with Hamiltonian Hq, controlled by an external parameter A (in 
our experiments, A is the longitudinal field B z ). We assume Hq to have gapped phases around 
a critical point A c , and without loss of generality we write a perturbed system Hamiltonian 
Hi — Hq + eV, where V is an arbitrary Hermitian operator (to be defined later) and e is 
the strength of the perturbation. Taking the ground state 1 ( A) ) of H as the initial state, 



the time dependent Loschmidt echo 27j takes the form 

L(t) = \£(t)\ 2 = | (0(A) | e iHlt e- iHot |0(A)> | 2 . (6) 

Notice that the evolution under Hq gives a physically irrelevant phase, which we keep for 
convenience of notation. The correspondence of the quantum critical points a QPT and 
the minima of the Loschmidt echo for long times has been shown for many systems 0, . 
However, the dynamical behavior for short times depends on the symmetries of the phases 
around the critical point and those of the perturbation operator. For instance, a monotonic 
increase of the decay rate with a singularity in its first derivative has been observed for some 



systems with local perturbations 28J. On the other hand, in the experimental results shown 
in the next section we observe that, for a fixed short time, the Loschmidt echo approaches a 
minima in the vicinities of the critical points. In this section we are concerned with providing 
a theoretical framework to these experimental observations. For this, we will analize the 
oschmidt echo for short times using a perturbative approach (similar to the one of Ref. 



281]). and particularize to the universality of the system we simulate in the experiments. 
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For small perturbations e we expand the echo amplitude 



e(t)~t(t)\ 



e=0 



+ 



d£(t) 



de 



d 2 £(t) 



£ = 



de 2 



£=0 



The first term is 



'(*)Uo = <0(A)| 



D iHot — iHot 



0(A)) = 1. 



(7) 



(8) 



For the second and third terms, we need to compute derivatives of the perturbed evolution 
operator. We can do this by expanding into infinite series and re-summing after computing 
the expectation value of the operators in the ground state. After some algebra (see appendix 
[£}, we find 

d£{t) 

£=0 



de 
d 2 £(t) 



-it)V c 



oo 



(9) 



de 2 



N-l 



£=0 



2j> 



a\ 2 x 



a=0 

-i(E a -Eo)t _ 1 + it ( Ea _ £ Q ) 



(B.-Btf • (10) 

where a indexes the N eigenstates of Hq with energy E a , Eq is the ground state energy, and 
V 0a = (a(A)| V |0(A)). The second order term of Eq. fflOl) resembles the so called fidelity 
susceptibility 29| and the quantum geometric tensor [30|] that have been shown to display 
singular behavior and scaling near a critical point. Indeed, if we take the Fourier transform 
of |L(£)| 2 , we obtain the fidelity susceptibility 29| for low frequencies. Higher frequency 
components appear that are related to the extra terms in the local density of states that 
generalizes the ground state fidelity {(J. 

Our final perturbative expression for the Loschmidt echo is then 

N-l 



L(t) 



a=l 



1 — cos(-E a — E )t 
{E a -E ) 2 



(11) 



B. Landau-Zener QPT toy model 



When the main contribution to the sum in Eq. (Ilip is given by the first excited state, 
we can approximate 



L(t) 



\Voi\ 2 
A 2 



e 2 (l 



cos At), 



(12) 



where A = E\ — Eq is the gap that has a minimum at the critical point, and we have assumed 
that there are no degeneracies. For degenerate systems like our experimental one, we just 



have to replace | Voi | 2 by a sum over the degenerate subspace of the transition elements 
squared. In a typical second order QPT, A ~ |A — \ c \~ zv , where v is the correlation length 
critical exponent and z is the dynamical critical exponent. As described in Sec. II, for a 
finite system the gap does not close but reaches a minimum A min that goes to zero with 
the size of the system N. Thus, non-analyticities occur only in the thermodynamic limit 
iV -> oo. 

Eq. ffl2|) suggests that whenever the ground and first excited states are the most relevant 
for a particular system dynamics, we can study the qualitative features of a QPT with a 
two-level toy system under both transversal and longitudinal fields, 

H LZ = A min a x + s(X)\X\ zu a z , (13) 

where s(A) is the sign function. Furthermore, this toy model - which represents Eq. f TT2"j) 
exactly up to 0(e 2 ) - resembles the approximations we use to model a QPT with our 
NMR quantum simulator - see Fig. (T5]) for a comparison between exact results and this 
approximation. 

From the spectra of our numerical simulations (see Fig. [2]), we see that our experiments 
are best described by zv — 1. In this case, Eq. ffl3|) is the well known Landau-Zener model 



3l| that has been used successfully to predict the scaling laws for the creation of topological 



defects when a system is quenched at finite speed through a critical point [32 ] . For this 
Landau-Zener model, 



A = 2 v /A2 + AL ri (14) 
A 2 

IT/ 1 2 _ mm ( 1 z.\ 

\ v oi\ ~ X2 — TT2' ^ 

min ' 

which means it has a "critical point" at A = 0. Expanding Eq. (fl2j) for short times, and 
replacing with Eq. ffl5|) . 

^M-sHW (16) 

Since the decay rate of L(t) (proportional to | Voi | 2 ) has its maximum at A = A c = 0, then 
we conclude that the decay of the echo is strongest at the critical point - or, conversely, that 
for a fixed time t the echo has a minimum at the critical point. 

In order to discuss possible universal scaling properties of the Loschmidt echo, our gener- 
alization in Eq. (fTB"]) from a Landau-Zener model attempts to incorporate a gap that closes 
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with an arbitrary power zv ^ 1. In this general case the short time decay is still given by 
Eq. f fT6l) . with a decay rate £ 2 |Voi| 2 . By choosing V = a z independent of A and e, and taking 
A m j n = 1/N for demonstrative purposes, we find that near the critical point, 

I V 01 1 2 ~ . (17) 

This suggests that the decay rate of the Loschmidt echo might show scaling with universal 
exponents. Such scaling has been proven for the ground state fidelity and the quantum 



geometric tensor 



30] . In principle, our experimental technique could be used to test univer- 



sality and scaling properties of the system. However, our experiments are currently limited 
to the case zv — 1 and relatively small sizes that prevent us from exploring these properties. 



IV. NMR IMPLEMENTATION 
A. Overview of the experiment 

Our goal is to measure the Loschmidt echo in the antiferromagnetic spin chain described 
by Hamiltonian |1| as a parameter (B z ) is varied, and from this infer the critical points 
of the system. Step by step, the experiment can be summarized as follows: Starting from 
the thermal equilibrium state, we prepared a pseudo-pure state (gate sequences for this are 
shown in Figures [3] and [7j) . For each value B z we transform the pseudo-pure state from the 
computational basis to the ground state of the Hamiltonian ([T|) (that depends on B z ) using 
a unitary Uq. We evolved the system forward in time with the Hamiltonian ([T|) at field B z , 
and then backwards with a perturbed field B z + e. After transforming the state back to 
the computational basis using Uq, we encode L as the diagonal element that is indicated 
by the initial computational basis in the density matrix. Exploiting another operation D to 
eliminate the non-diagonal elements of the density matrix, we can obtain the locations of 
the minima of L using a selective readout pulse and observing the intensity of a spectrum 
of a single qubit. We perform the experiment in chains of three and four spins. The results 
are shown in Figures 151 and [TUj respectively. 

We simplify the implementation of the experiment with a number of approximations 
summarized here and described in detail in the following sections. At each value of B z we 
prepare a very good approximation of the ground state, with fidelity higher than 98% (we 
elaborate on this point in the conclusions). The approximated ground state is obtained with 

9 



perturbation theory around the crossover point of zero transverse field and does not require 
knowledge of the criticality of the system with non-zero transverse field. We split the range 
of the field B z in intervals (three and four for the odd and even spin chains respectively) 
and use a different pulse sequence for each interval. The forward-backward evolution is 
compressed into a single step using a first order Trotter expansion with 98% accuracy. The 
quantum networks for the odd and even chain experiments are shown in Figures H] and [HJ 



B. Efficient detection of critical points using the Loschmidt echo 

In order to measure the Loschmidt echo we first prepare the ground state \tp(B z , B x )) 
of H(B Z , B x ), which remains very close to Eqs. ([2]) and (j4j), except in the vicinity of the 
critical points. Then, we evolve it forward under H for a period of time t, and next evolve it 
backwards under H + eV for t, where eV is the fixed perturbation introduced for detection 
with |e| <C 1. Here, B z will be our control parameter (A in the previous section), and we 
choose the perturbation as V = — YliLi a l- This choice of a global perturbation simplifies 
our experiments, although more general choices like local perturbations lead to the same 



results but with a reduced signal 



2s 



In order to detect the critical point of the transition 
we fix the evolution time t = t and the transversal field B x , and measure L as a function of 
B z f], S|- As shown in the previous section, the critical points will be marked by the minima 
of 

L = L(B z )\ t=r = mB z ,B x )\UlU\ij(B z ,B x ))\ 2 , (18) 

where U = e~ lrH and U p = e- tT ^ H+£V ) are the unperturbed and perturbed evolution opera- 
tors, respectively. We show some representative echoes in small chains in Figure [2](c-d). 

Measuring an overlap such as Eq. (fl8|) in general might require full state tomography 
techniques. Because of its particular form, we can also couple the system to a probe qubit 
in such a way that L is encoded in the off-diagonal terms of the reduced density matrix of 
the probe ^Js)]. Here, we present a method to measure L directly in the system. We call U 
the unitary operation that prepares \tp g (B z , B x )) from an arbitrary computation basis state 
\s). This is not necessarily an efficient operation for all systems - indeed, finding the ground 



state of arbitrary Hamiltonians might be an NP-hard problem 



However, theoretical 



results suggest that 



any initial state with a large overlap with the ground state is sufficient 



to detect criticality 34|. For instance, in our experiments we do not prepare the true ground 
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state of the system, but actually a state that approximates it very well. We will discuss this 
and other alternatives to the preparation of the ground state in the last section. 
Through rewriting Eq. (ITS)) as 

L=\{8\UtU}UU \ S )\ 2 , (19) 

we find that L can be obtained by projecting 

|tf> = f^f/f/ols) (20) 

onto state |s), i.e. L is equal to the element \s)(s\ of the density matrix p = \^f)(^\. Without 
loss of generality, we chose \s) = |00...0), the state with all qubits in computational basis 
state |0). After the final evolution Uq, we eliminate the non-diagonal elements by gradient 
pulses or dephasing processes [35|, |36(. Then, through a read-out pulse, e.g. tt/2, applied to 
an arbitrary qubit, we obtain the signals marked by the states of other qubits. 

We are only concerned with the signal marked by the state in which all other qubits are 
in state |0). Because in NMR we observe differences in populations, the amplitude of this 
signal A is proportional to (L — p nn ) < L, with 1. The locations of the minima of A are 
the same as those of L, with their values each decreased by an additional p nn . This allows 
us to detect the critical points through A by measuring only one qubit of the system. 

C. Odd N case 

We first demonstrate the detection of critical points of a QPT in an odd spin system with 
N = 3. We prepared an initial state that approximates the ground state of the Hamiltonian 
for each value of B z . Using our notation for the ground states of H for B x = (|^) = 1 000) , 
|010), 1 101) , and 1 1 1 1) , for k = 1..4 respectively), the ground state near B c = ±2 can be 
approximated as 

mB z ,B x )) = \rjcos<p-\r n )sux<p (21) 



with 



tany? 



[2 - \B Z \) + ^{2-\B z \f + Bl\ /B x , (22) 



where m — 1, n — 2 or m — 4, n = 3, corresponding to B c = — 2 or 2, respectively. In 
the vicinity of B c = 0, the gap between the lowest energy levels is so small that the ground 
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state can be well approximated by \ip 2 ), (IV^) — / or IV'I); corresponding to B z < 0, 
£>z = 0, or B z > 0, respectively. 

For the experimental implementation, we used 13 C labelled trichloroethylene (TCE), 



dissolved in d-chloroform as the sample 37| • Data were taken with a Bruker DRX 700 MHz 
spectrometer. We denote the X H nuclear spin as qubit 2 (H2), the 13 C directly connected to 
1 H is denoted as qubit 1 (CI), and the other 13 C as qubit 3 (C3). The difference of frequency 
between CI and C3 is about 1249.2 Hz, and the coupling constants are J13 = 103.1Hz, 
J12 = 200.9Hz, and J 23 = 9.16Hz. The spin-selective excitation for CI or C3 is realized by a 



GRAPE pulse 



38 



The J-coupling evolution e l< t >a -^ cr ~ between qubits I and k is implemented 



by a standard refocusing pulse sequence [39j. The effective pure state 1 000) is prepared 
by spatial averaging [40( from the thermal equilibrium state p eq = Jh&z + lc{p\ + cr z ), 
by approximating the system as a weakly-coupling system, where 7# and 7c denote the 
gyromagnetic ratios of proton and carbon, respectively. The gate sequence for the pseudo- 
pure state preparation is shown as Figure [3j 

In order to measure the echo we split the B z axis in intervals near B c = —2, 0, and 2. In 
particular, we use different quantum networks for B z G [—3, —1], (—1, 1), and [1, 3], show in 
Figures H] (a-c) respectively. The operations for preparing Uq and Uq are indicated by the 
dashed rectangles and D denotes the operation to eliminate the non-diagonal elements of 
the density matrix. Figure [5] shows the corresponding gate sequences. The evolution time 
is chosen as r = tt, and e = 0.2 or 0.125. The echo evolution U^U can be approximated 
by UpJ ~ e~' LT£ ^ +CT 'z + ' T ^ with fidelity larger than 98%. We optimize the gate sequence 



4l| 



CNOT 21 - e~ iT£a * - CNOT 21 as e ~ iT£a ' a ", and CNOT 23 - e- iT£a * - CNOT 23 as t~ iT£ ^ 
to obtain figure 0(b) from Figure H](b). The amplitudes of signals are obtained by measuring 
on H2, with experimental results shown in Figure Experimental data are marked by " x" 
and "+" for e = 0.2 and e = 0.125, respectively. The corresponding theoretical results are 
indicated by the light and dark curves. The critical points are correctly indicated by the 
minima of the amplitudes. 

D. Even N case 

We illustrate the detection of QPT critical points in an even spin chain with N = 4. Here 
we use the notation for the ground states for B x = 0, \tp e k ) = |0000), (|0100) + |0010>)/ v / 2, 
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(|0101) + |1010))/v^, (|1101) + |1011))/v^, and |1111), for k = 1.. 5 respectively. Depending 
on the value of B z we prepare an approximation to the ground state 

t/j(B z ,B x ) = \rjcos<p-\r n )sm<p. (23) 

For B, near B r = ±2 



tan<p = [(2 - \B Z \) + y/(2 - \B z \f + 2BI }/{V2B x ) (24) 

with m — 1, n — 2, or m — 5, n — 4, corresponding to I? c = —2 or 2 respectively. For S 2 
near B c = ±1 we use 

tan^= [(1 -[£,[) + 7(1 -\B z \f + Bl)/B x (25) 

with m = 2, n = 3 or m = 4, n = 3, corresponding to 5 C = — 1 or 1, respectively. 

For implementation, we choose the four carbons in crotonic acid 42j dissolved in d6- 
acetone as the four qubits by decoupling the protons. Data were taken with a Bruker 
DRX 700 MHz spectrometer. The chemical shifts for the four carbons z/i_4 are —2965.75, 
-25501.9, -21583.9 and -29431.5 Hz. The J-couplings are J 12 = 41.6, J 23 = 69.7, J 34 = 
72.0, J13 = 1.5, J u = 7.0, and J 24 = 1.2 Hz. 

We prepare the pseudo-pure state by spatial averaging through improving the scheme 



found in [43j. Our technique can be illustrated by transforming the thermal equilibrium 
state of a four qubit system Ylt=i a l ^° 

^^(1/2 + ^+^/8, (26) 

where 1 denotes the unit operator and Ylt=i a \ can be transformed to an effective pure state 
in the three-qubit system. This method generalizes to an iV-qubit system in a recursive 
manner. After some simplifications [y], the complete gate sequence to generate |0000) 
is shown as Figure [71 where the state specific swap gate requires two J— couplings with 



evolution time 1/(2^) |45j. In the ideal case the strength of the single peak obtained 
through a 7r/2 read out pulse selective for one spin is equal to that of the same peak in a 
spectrum of the thermal state, where eight peaks with equal strength appear. 

The ground states are prepared from Eqs. ([23]- [25]) . As before, we split the B z axis 
in intervals around the critical points of zero transverse field. The networks to measure 
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the echo for B z £ [—3,-1.44], and (—1.44,0] are shown in Fig. [8] (a-b). From these one 
can obtain the networks for the intervals B z £ [1-44, 3] and (0, 1.44) simply by adding 
NOT gates to all qubits at the end of the corresponding networks for implementing Uq. 



Through compiling the pulse sequence 11|, we obtain the gate sequences shown as Figure El 
where UpJ ps e~ lT£( - cr z +cr z +cr z +cr ^ with fidelity larger than 98% and the two SWAP gates are 
cancelled because they commute with e _rre(c ^ +<J ^. Experimental results are shown in Figure 
[TUl with t = 7r/2. The measured amplitudes are marked by " x" and "+" for e = 0.5 and 0.4, 
respectively. The solid curves shows the corresponding theoretical results. Again, the critical 
points are correctly indicated by the minima of the amplitudes, so the experiment results are 
in good agreement with theoretical expectations. The observed errors could be explained 
by imperfections in the implementation of the radio frequency pulses, inhomogeneities of 
magnetic fields and decoherence. 



V. DISCUSSION AND CONCLUSIONS 



We performed experimentally quantum simulations of the second order quantum phase 
transitions in finite systems. In particular, we showed the QPTs and found the critical points 
of three- and four-spin Ising chains, representative of odd and even spin chains, respectively. 
The critical points are indicated by the minima of the Loschmidt echo. We showed that this 
echo can be realized by inducing the perturbation with an external field, and the positions 
of its minima (related to the critical points) can be obtained by measuring only an arbitrary 
qubit of the system. In the weakly and fully resolved coupling systems, the resonance lines 
can be assigned, and the line marked by the other qubits in |0) can be identified. However, in 
large size systems where the requirement of fully resolved couplings is not practical, or in the 
strongly coupling systems, e.g., liquid crystal or solid NMR systems, where the assignment 
of resonance lines are not possible, one cannot identify the marked line. For these cases, 
our method can be generalized by measuring the global polarizations of the whole system 
by a collective n/2 pulse (or N pulses selective for each qubit), replacing the readout pulse 
applied to one qubit. In the vicinities of the critical points, the loss of the polarization 
due to the decoherence process (e.g., gradient pulse or dephasing process) approaches the 
maxima. Hence the critical points will be indicated by the minima of the amplitude of the 
total signals of all qubits. Furthermore, this has the advantage that a global measurement 
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is scalable with the size of the system. 

Our method improves the previous one that required a probe qubit for both the per- 
turbation and the measurement j^j. We believe this advantage gives our method better 
scalability with the size of the system. In particular, the perturbations created by the 
probe qubit method are limited by the probe-system coupling strength, and, furthermore, 
can become weaker than the noise in large systems when they do not couple the probe to 
a macroscopic number of normal modes in the system. Separating the perturbation and 
measurement also gives finer control over the whole experiment. 

On the issue of scalability, a very important point in the algorithm is the preparation of 
the initial state. From a theory point of view, most of the studies of the Loschmidt echo 
have used ground states as initial states only because of simplicity. However, preparing the 



ground state of an arbitrary Hamiltonian is an NP-hard problem |33|]. Furthermore, it would 
be redundant, since it is most likely that knowing the exact ground states is equivalent to 
knowing everything about the system - including the information about criticality that one 
wants to obtain from the echo experiments. Nonetheless, there is evidence that the initial 
state need not be the exact ground state, but any state with a sizeable overlap with the 
ground state. For instance, analytical studies show that thermal states at temperatures 
at or below the energy scales of the system can be used effectively to detect the quantum 



phase transition 



34|, where the number of the spins can be up to 200. However, in some 



systems (like our liquid NMR experiments) preparing a thermal state is not particularly 
easier than other -perhaps more useful- states. For instance, in our experiments we prepared 
a good approximation to the ground state that we obtained from a simple perturbation 
theory around the crossover point of zero transverse field. This method suggests that other 
approximations, such as mean field or numerical classical algorithms, could work well to 
detect criticality. 

While the problem of finding strict minimum requirements for the initial state of the 
algorithm is clearly in need of more research, we feel that it is reasonable to argue that 
initialization of the algorithm is scalable: it only requires finding among many possibilities 
one that can be prepared efficiently in a quantum computer. It would be interesting to study 



the effect of more efficiently-prepared ground states [33( or to investigate if state-independent 



indicators - such as the operator fidelity susceptibility proposed in Ref. 



461 ] - could get rid 



of the initial state issue altogheter. Finally, we would like to mention that other possible 
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extensions of our experimental methods are using the Loschmidth echo to measure QPTs in 
gapless systems and also for measuring thermal phase transitions [29|, |48j . 



VI. ACKNOWLEDGEMENTS 

We thank D. Suter, J. P. Paz, C. Batista, G. Ortiz, T.-C. Wei, M. Zwolak, and H.T. Quan 
for helpful discussions. F.M.C. acknowledges support from DOE-LDRD programs. J.K.G. 
acknowledges support from the NSF-REU and Los Alamos Summer School programs. 



APPENDIX A: PERTURB ATIVE EXPANSION OF THE LOSCHMIDT ECHO 



We start from the expansion of Eq. (J7j) , 

d£{t) 



de 



e + 



d 2 i{t) 



e=0 



de 2 



(Al) 



e=0 



where 



l{t) = (0(A) | e ^ t e -i(H 0+ eV)t | Q(A) ^ (A2) 

with 1 0(A) ) the ground state of H , and we keep the harmless e tHot operator because it will 
simplify the results. The first term of the expansion can be simply evaluated as in Eq. ([8]), 



^(*)Lo = (°( A )I e iHot e- tHot |0(A)) = 1 



(A3) 



For the first and second order terms we must compute derivatives of the evolution operator. 
We can do this by expanding the exponential into an infinite series sum, 
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Computing now the expectation value 
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(A5) 



where E is the ground state energy and Voo = (0(A) | V 1 0(A)) . 

For the second order term we continue derivating Eq. (IA4I) before the evaluation at e — 0, 
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By taking the expectation value on the ground state we now obtain 
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replacing now k' = k — 1 and m! = m + k in the first and second sums inside the brackets 
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We can simplify the term inside the brackets by counting the number of times the terms 
with \k — m\ = 0, \k — m\ = 1, and so on are repeated. The final expression is then 
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We can make further progress by inserting identities ^2^=o \ a ) ( a \i with \a) the basis of 
eigenstates of H (we assume a finite Hilbert space a = 0, .., N — 1), 
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where Vo a = (a\ V 1 0(A)) . We can do the sum over k first, 
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(notice that the term with a = is finite), followed by the sum over n, 
d 2 £(t) 



de 2 



N-l 

e=0 a= Q 



e -i(E a -E )t _ 1 + ^(^ _ ^ 



'Oct I 



7V-1 



\v 0Q \H 2 -2^2\v Qa \ 



a=l 



1 _ e -»(g a -g )t _ it ( Ea _ ff Q ) 

(E Q - E ) 2 



(A12) 



Now we need to put the results of Eqs. (jA~3ft . flA5j) . and flA12j> into Eq. (jATJ), 
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Using that Voo is real and keeping the term with lower order in e, we obtain the expression 
for the Loschmidt echo: 

L(t) = \l(t)f * 1 - 2e 2 £ \V o y-^-f o)t (A14) 
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FIG. 1: Phase diagram of the antiferromagnetic Ising chain with transverse and longitudinal fields, 
B x and B z respectively, in the thermodynamic limit of infinite chain size 
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241 ] . The coupling 



strength is chosen as the unit for B x and B z , [see Eq. ([I]) in text]. In the shadowed region 
inside the circle the ground state is (doubly degenerate) antiferromagnetic (AF), and in the clear 
region outside it the ground state is paramagnetic (PM). The transition line between both phases 
is a second order critical line, while the points at B x = are first order transitions. The phase 
diagram corresponds to that of a two dimensional classic Ising model with field equal to B z and 
effective temperature proportional to B x . The dashed line shows qualitatively the region we explore 
experimentally, where the critical points in the thermodynamical limit are close to B z = ±2. For 
the finite systems used in our experiments, we need to consider boundary effects which show up 
like extra sub-phases inside the AF phase. For odd N, a new critical point appears at B z = 0, 
while for even N two extra critical points appear at B z = ±1. 
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FIG. 2: Phase diagrams without transverse field (a-b) and Loschmidt echo (dimensionless) with 
small transverse field (c-d) for the Ising-chains with odd and even spins, shown in the left and right 
columns, respectively. The dark and light curves in figures (a-b) represent the two lowest energy 
levels, by setting the coupling strength and Ti to unity. The phases and energy levels are listed in 
Eqs. ([2]E]). The crossover points are B c = ±2, in the odd spin system, and B c = ±2, ±1 in the 
even spin system. The minima of the Loschomidt echo in panels (c) and (d) indicated the critical 
points. Without loss of generality, we choose N = 7 and 8 to illustrate the odd and even cases, 
where e = 0.1, r = it, and B x = 0.1, for calculating L. In figures (c-d) the light thick curves show 
the numerical results from Eq. (J6j), while the dark thin curves show the approximate analytical 
results from Eq. (1121) . 
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FIG. 3: Gate sequence to prepare the effective pure state |000) by spatial averaging from thermal 
equilibrium state of TCE, where cos a = 2^c/lH- Here "yn and jc denote the gyromagnetic 
ratios of proton and carbon, respectively. The single qubit gates are implemented through radio 
frequency pulses denoted by the rectangles. The rotation angles and directions are shown inside 
and above the rectangles. The bold vertical lines denotes the gradient pulses along z axis. The 
two filled circles connected by a line denote the J— coupling evolution e _ * <? " Jz<Jz between qubits I 
and k, where <p is shown next to the line. 
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FIG. 4: Quantum networks for measuring critical points in intervals B z £ [—3,-1], (—1,1), and 
[1, 3] in the three qubit system, shown as figures (a-c) respectively. R = e l<firTy where ip is given by 
Eq. (f22j) . and i?o = 1 (unit operator), e" 17 "/ 4 or e l7r<Jy / 2 for B z = -0.5, 0, or 0.5, respectively. Uq 
and C/g are indicated by the dashed rectangles, and UpU ~ e _ * T£( - <T ^ +<J ^ +<J ^. (J) and the black dot 
connected by a line denote a controlled NOT gate, and N denotes a NOT gate. D denotes the 
operation to eliminate the non-diagonal elements of the density matrix. The last operation in each 
figure denotes the measurement, which can be applied to an arbitrary qubit of the system. 
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FIG. 5: Gate sequences (a-c) to implement Figures E] (a-c), respectively. A 16-step average over a 
random delay, denoted by d, between and 10 ms, dephases the residual zero-quantum coherence. 
The last ir/2 pulse is the readout pulse, which can be applied to an arbitrary qubit of the system. 
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FIG. 6: Experimental results in the three qubit QPT system, where r = tt. The four phases 
\tp%) with k = 1, 4 are represented as = |000), 1 010) , 1 101) , and 1 1 1 1) , respectively. The 
experimentally measured amplitudes of the signals are marked by "x" and " + " for e = 0.2 and 
0.125, respectively. The minima of the amplitudes indicate the critical points. The theoretical 
results are shown as the light and dark curves. The experimental results show a good agreement 
with theory. 
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FIG. 7: Gate sequence to prepare the effective pure state 1 0000) by spatial averaging from the 
thermal equilibrium state of the four carbons in crotonic acid, where cos a = 1/8, and cos/3 = 1/4. 
The filled rectangles in pairs connected by a line denote a state specific swap gate between qubits 
I and k, i.e., it transforms a l z to cr^, or to a\. 
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FIG. 8: Quantum network for measuring critical points in intervals B z € [—3, —1.44], and (—1.44, 0] 
in four qubit system, shown as figures (a-b) respectively. H denotes the Hadamard transform gate, 
and u}ll « e " ire ^ +fT ' +<T " +fT ^. R 2 = e itfl2 ^ and R x = e iipia y , where <p 2 and tpx are chosen as 
Eqs. (|24H25p . The rectangle and the dot connected by a line denote a controlled operation that is 
shown inside the rectangle. The filled rectangles in pairs connected by a line denote a SWAP gate. 
The networks for intervals B z £ [1.44, 3] and (0, 1.44) can be obtained by adding NOT gates to all 
qubits at the end of the networks for implementing Uq in figures (a-b), respectively. 
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FIG. 9: Quantum gate sequences (a-b) to implement Figures [8] (a-b), respectively. Through re- 
placing e by — e in figure (a-b) one can obtain the gate sequences for the intervals B z G [1.44,3] 
and (0,1.44), respectively. 
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FIG. 10: Experimental results in the four qubit QPT system, where r = tt/2. The five phases |^|) 
with k = 1, 5 are represented as \ip e k ) = |0000), (|0100) + |0010))/\/2, (|0101) + |1010»/\/2 5 
(|1101)+|1011))/v / 2, and 1 1111), respectively. The experimentally measured amplitudes are marked 
by " x" and "+" for e = 0.5 and 0.4, respectively. The minima of the amplitudes indicate the critical 
points. The theoretical results are shown as the light and dark curves, in good agreement with the 
experimental results. 
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